Quantum information processing based on 31 P nuclear spin qubits 
in a quasi-one-dimensional 28 Si nanowire 
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We suggest a new method of quantum information processing based on the precise placing of 31 P 
atoms in a quasi-one-dimensional 28 Si nanowire using isotope engineering and neutron-transmutation 
doping of the grown structures. In our structure, interqubit entanglement is based on the indirect 
interaction of 31 P nuclear spins with electrons localized in a nanowire. This allows one to control the 
coupling between distant qubits and between qubits separated by non-qubit neighboring nodes. The 
suggested method enables one to fabricate structures using present-day nanolithography. Numerical 
estimates show the feasibility of the proposed device and method of operation. 



I. INTRODUCTION 

There is great scientific and commercial interest in the 
development of quantum computation (QC) and the cre- 
ation of computational devices based on the principles 
of quantum logic. Several different schemes for QC have 
been proposed to date (see, for example, Ref. 1 and ref- 
erences therein). One of the exciting avenues, potentially 
compatible with the vast fabrication capabilities of mod- 
ern semiconductor technology, relies on the encoding of 
information in the electron or nuclear spins present in 
semiconductor nanostructures, leading to a spin-based 
semiconductor quantum computer [2-6] . 

The most developed model of the nuclear spin quantum 
computer is the Kane suggestion [4] to use a precisely 
located array of phosphorous donors introduced into Si. 
In this proposal, the nuclear spin l J% of 31 P is used as a 
qubit, while a donor electron together with an overlying 
gate (A-gate) separated from the donor by a SiC>2 or 
Sio.85Geo.15 barrier, provides single-qubit operation using 
an external magnetic field and pulses of radio-frequency 
radiation. The interqubit coupling is determined by the 
overlap of the electron wave functions and is controlled 
by metallic gates (J-gates) midway between the A-gates. 
The overlap of wave functions of localized electrons in Si 
drops very rapidly with distance r, exp(— 2r/ae), where 
Ob is the radius of localization (for P in Si, ae ~ 2.5 nm), 
therefore the interqubit distance r between P atoms must 
be small (less than 20 nm) to allow overlap. 

Experimental realization of the suggested model 
presents a number of difficulties. We focus here on two 
problems: (1) placing single P donors into the Si sub- 
strate at a precise depth underneath the barrier and (2) 
the necessity to increase significantly the interqubit dis- 
tance r to have room enough to arrange the metallic 
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gates. This means that mechanisms other than the direct 
overlap of the electron wave functions have to be chosen 
for the coupling of adjacent nuclear spin qubits. 

To solve these problems, we suggest the novel technol- 
ogy based on epitaxial growth of Si and SiGe layers from 
isotopically engineered Si and Ge sources followed by 
neutron-transmutation doping of the grown structures. 
We also describe the mechanism of indirect interqubit 
coupling based on the arrangement of qubits in a meso- 
scopic quasi-one-dimensional wire. This mechanism al- 
lows one to control the coupling between qubits sepa- 
rated by large distances r (200 nm or even more), which 
permits the fabrication of metallic gates by means of 
the modern lithography. Moreover, the suggested mech- 
anism of indirect coupling allows entanglement passing 
over non-qubit nodes in an array of qubits. We also 
present the numerical estimates which justify the feasi- 
bility of the proposed device and method of operation. 



II. DEVICE FABRICATION 

A. Precise placement of P atoms into Si 

Two methods had prevoiusly been suggested for the 
precise placing of P atoms into Si (for a review, see 
Ref. 7 and references therein). In [8] and [9], a "bottom- 
up" method is described for incorporation of phosphorus- 
bearing molecules PH3 on a preliminary H passivated 
Si (001) surface during the molecular-beam-epitaxial 
growth followed by the decomposition of PH 3 at sig- 
nificantly increased temperature. The alternative "top- 
down" method is based on incorporation of dopant atoms 
under the surface of the grown structure using ion im- 
plantation, followed by the annealing of radiation damage 
[10-12] at increased temperatures. In these methods, the 
depth distribution of the incorporated donor atoms can- 
not be controlled to the necessary accuracy. For example, 
in the "bottom- up" method, incorporation of impurities 
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FIG. 1: Schematics of the proposed device. After NTD, 31 P donors appear only inside the 30 Si-spots and underlying 74 Ge-strips 
will be heavily doped with 75 As donors. All sizes are shown in nm. 



of different chemical nature is disadvantageous both for 
the quality of the growing structure and for the sharpness 
of the vertical distribution of impurities because of the 
" floating- up" effect in the process of subsequent growth. 
In the "top-down" method, the distribution of impurities 
in " as-implanted" samples is determined by the energy of 
ions, and initially is not sharp. Moreover, during the an- 
nealing of radiation damage, impurities are redistributed 
towards the Si/SiGe interface which acts as a getter. As 
a result, the vertical distribution of the implanted impu- 
rities is rather extended. Meanwhile, small fluctuations 
in the vertical position of P atoms under the gate is very 
important to minimize the A-gate voltage error rate [13] 
because otherwise, each qubit would need its own set of 
applied voltages and each interacting pair of qubits would 
need its unique pulse time [14]. 

We suggest here a novel technology of the precise plac- 
ing of P atoms into Si layer. The key point is the growth 
of the central Si and barrier Sio.85Geo.15 layers from dif- 
ferent isotopes: the Sio.85Geo.15 layers from isotopes 28 Si 
and 72 Ge and the central Si layer from isotope 28 Si with 
30 Si spots introduced by means of the nano-lithography 
(Fig. 1). The formation of quasi-one-dimensional Si wires 
will be achieved in a subsequent operation by the etching 
of Si layer between wires and the filling of the resulting 
gaps by the Sio.85Geo.15 barrier composed from isotopes 
28 Si and 72 Ge. Because different isotopes of Si and Ge 
are chemically identical, this technology guarantees the 
high quality of the grown structures. 

After preparation, these structures will be irradiated 
with a neutron flux in a nuclear reactor followed by the 
fast annealing of radiation damage. The behaviour of 
different isotopes is different. After capture of a slow 
(thermal) neutron, a given isotope shifts to the isotope 
with mass number larger by one. If the isotope thus 
obtained is stable, this nuclear reaction does not entail 
doping. However, if the obtained isotope is unstable, it 
transmutes after half-life time r to a nucleus of another 
element with atomic number larger by one in the case 



of (3~— decay. This method of doping is called NTD 
- neutron-transmutation doping [15]. In the case of Si, 
NTD is based on the transmutation of the isotope 30 Si: 

3 °Si + In = -+0-(t = 2.62/t) -> HP. 

In the isotopically engineered structure, 31 P donor 
atoms will be produced only within 30 Si spots, because 
the isotopes 28 Si and 72 Ge shift to the stable isotopes 
29 Si and 73 Ge, respectively. Therefore, in the suggested 
method, the processes of the structure growth and dop- 
ing arc completely separated. The idea to fabricate a 
basic element for a nuclear spin quantum computer us- 
ing the isotope engineering of Si and SiGc nanostructures 
was proposed earlier in [16, 17]. The precise placing of 
P atoms into a Si matrix by means of the NTD method 
was suggested in [18]. 

Let us estimate the feasibility of the proposed method. 
We consider a 30 Si spot of area S = 30 x 30 nm 
w 10~ n cm 2 , thickness d =10 nm with a distance of 
200 nm between spots (Fig. 1). The "buried" distance to 
the structure surface is halved, 100 nm, which is needed 
to protect against cross-links and to ensure the influ- 
ence of the A-gate voltage on the corresponding under- 
lying qubit only. In the proposed method of incorpora- 
tion of P into Si, the vertical accuracy of the location 
of P donors is determined by the thickness of the 30 Si 
spot (10 nm) with respect to the distance to the A-gate 
(100 nm). Indeed, the irradiation of samples by ther- 
mal neutrons occurs at room temperatures where the in- 
troduced P atoms are immobile. The irradiation is fol- 
lowed by the annealing of radiation damage at higher 
temperatures (700°C). During the short annealing time 
(1 hour), impurities cannot diffuse far from the transmu- 
tation place; the diffusion length docs not exceed 1-2 nm. 
As a result, the proposed method will provide an almost 
equal burying depth of P atoms with a controlled accu- 
racy (about 10% in our example). The time needed for 
the irradiation of the grown structures in a nuclear reac- 
tor is estimated as follows. The number of transmutation 
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FIG. 2: Probability to find spots with 0, 1, 2 and 3 impurity 
atoms, introduced by NTD, as a function of the dimensionless 
irradiation dose $*. 

events is iV( 31 P) = 7V( 30 Si)cr 30 <I>, where N is the num- 
ber of atoms, a is the cross-section of the thermal neu- 
tron capture for given isotope (a 30 « 0.11 • 10~ 24 cm 2 ), 
$ = ipt is the integral neutron flux (</? is the intensity 
of the thermal neutron flux and t is the time of irradi- 
ation). In a spot with volume V — 30 x 30 x 10 nm 
ps 10~ 17 cm 3 , there are N = N(Si)V « 5 • 10 5 Si atoms 
(JV(Si) = 5 • 10 22 cm" 3 ). If the enrichment of Si with 
isotope 30 Si is close to 100%, N( 30 Si) ps 5 • 10 5 . To 
achieve 7V( 31 P) = 1, the integral irradiation dose $ has 
to equal to 2 • 10 19 neutron/cm 2 . In some research nu- 
clear reactors, <p ^ 10 14 cm -2 s , which corresponds to 
the reasonable irradiation time t ~ 2 ■ 10 5 s w 56 hrs. 



B. Qubit certification 

An unavoidable peculiarity of the NTD is the casual 
character of the neutron capture. As a consequence, after 
NTD, some of 30 Si spots will contain no donor atom ("0- 
spot") and cannot serve therefore as qubits, while some 
will contain more than one donor atom. The probability 
P m to find "0-spot" , "1-spot" , "2-spot" , and so on (to = 
0, 1, 2, ...), is described by the Poisson distribution: 

P m = C% P m (l-p) D - m 

where is the binomial coefficient, p — A^( 30 Si)<7 30 <i 
is the probability for the neutron to be captured in a 
layer of thickness d, D = is the dimensionless dose 
of irradiation. 

These probabilities are plotted in Fig. 2 as a func- 
tion of the dimensionless parameter <I>* = pD. The best 
candidates for qubits are "1-spots" , while spots having 
more than one donor could also be considered as qubits 
after corresponding theoretical investigation. Only spots 



without donors ( "0-spots" ) are obviously non-qubits. In 
the NTD method, the maximal portion of "1-spots" is 
37% at $* = 1. In this case, the portion of non-qubit is 
also about 37%. If "2-spots" are also considered as pos- 
sible qubits, the optimal integral dose will correspond to 
$* = 2. In this case, the fraction of non-qubits decreases 
tow 14%. 

To determine the number of donors in each spot, we 
propose to use narrow source-drain (SD) channels fab- 
ricated beneath each donor position (Fig. 1). If the 
given spot contains one donor, it will form (together 
with the underlying SD channel and overlying A-gate) a 
flash-memory field-effect transistor (FET) with the qubit 
donor acting as a "floating gate". The donor electron 
is separated from the SD channel and A-gate by the 
Sio.85Geo.15 barriers of about 10-20 nm width and 100 
meV height [5] which are used for the electron confine- 
ment. However, a pulse of relatively strong voltage ap- 
plied between the A-gate and SD channel will tilt the 
barriers leading to electron transfer and donor ionization. 
It was shown in Ref. 19 that if the dopant-gate separa- 
tion distance is more than 8ob ~ 24 nm for P in Si (in 
our case this condition is satisfied because the distance to 
the gate is about 100 nm), the transfer of electron from 
donor takes place abruptly at a threshold voltage. An 
electric field of the positively-charged donor ion located 
only 10 nm from the narrow SD channel will shift the 
FET cutoff voltage at the value of about 10 mV, which 
is easy to observe [20] . If there are two or more donors 
in the given spot, the cutoff shift will be even larger. If 
there is no donor underneath the given gate, the shift will 
not be observed. 



C. Read-out operation 

We believe that the proposed SD channel can also be 
used for the read-out operation, i.e. for the detection of 
a single nuclear spin state. The direct control of a nu- 
clear spin state via nuclear magnetic resonance (NMR) 
measurements is a difficult problem. In Ref. 21, coherent 
control of the local nuclear spin was demostrated, based 
on pulsed NMR in a quantum Hall device. In Ref. 22, a 
self-contained semiconductor device is described that can 
control nuclear spins in a nano-scale region. Measure- 
ments of the electron spin state are much easier taking 
into account the possibility of a spin-to-charge conversion 
for electrons. In accordance with the Kane model [4], the 
state of the nuclear spin 31 P is mediated by the spin of 
donor electron via the hyperfine interaction. Therefore, 
the task is to determine the spin orientation of the cor- 
responding donor electron. The suggested method [6] is 
based on the fact that at low temperatures, a donor atom 
can capture the second electron with small ionization en- 
ergy, about 1 meV, which results in the appearance of 
a negatively charged donor (£> - -center). However, this 
process is possible only when the spin orientation of the 
second electron is opposite to that of the first electron. 
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The appearance of the charged donor in the vicinity of 
the narrow SD channel will affect the current [23, 24] and 
can therefore be detected. As a result, one can determine 
the spin orientation of two neighboring donor electrons 
if one applies a potential difference between the corre- 
sponding A-gates which will cause the electron to jump 
from one donor to another . If we choose the spin orienta- 
tion of the given donor as a reference, one can determine 
the spin state of the neighboring qubits on the right and 
left sides. 
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FIG. 3: Spatial distribution of electron wave functions for N 
electrons in a quantum wire of the length /. 



D. Source-drain channels fabrication 

The proposed SD channels have a twofold purpose: the 
determination of the number of donors within each spot 
(qubit certification) and the determination of the spin 
state of the donor electron (read-out operation). We sug- 
gest a method for the fabrication and the proper align- 
ment of the SD channels with respect to the position 
of 31 P qubits. This method is also based on isotope 
engineering of a Sii^Ge^ layer followed by NTD. Wc 
propose to make the underlying Sii^Ge^ layer from a 
composition close to pure Ge (say, Sio.iGeo.g) using iso- 
topes 28 Si and 72 Ge, followed by the fabrication of 30 
nm-width strips where 72 Ge is replaced by 74 Ge (Fig. 1). 
After NTD, these 74 Ge-strips will be doped by As donors 
through the following nuclear reaction: 

\\Ge + = ^Ge -> /T (r = 82 min) -> IjAs 

Irradiation of the structure with the thermal neutron 
integral dose <pt = 2 • I0 19 cm~ 2 (<f>* = I) needed for 
introducing on average one P donor in each 30 Si spot, 
will also lead to doping of 74 Ge-strips with As donors to 
a high level (ATa s ~ 4.5 • I0 17 cm -3 ) because of the rela- 
tively large <r 74 = 0.5 • 10~ 24 cm 2 . This concentration of 
As exceeds the critical value of the metal-insulator tran- 
sition for Ge:As [25]. Therefore, NTD-introduced narrow 
channels will have a metallic-like conductivity and remain 
conductive down to T — ► 0. This is important because 
nano-FET will operate at low temperatures when donor 
electrons in 31 P-qubits are localized on their donors. For 
the suggested geometry of the SD channel, with the thick- 
ness of the Sio.iGeo.g layer of about 10 nm, width of the 
74 Ge-strips of 30 nm and the length about 1 /im, the 
channel resistance is about 1 MCI, which is suitable. 

In the proposed method of device fabrication, the 
proper alignment of the 74 Ge-strips with respect to the 
overlying 31 P qubits is provided by the high accuracy 
of the electron beam-assisted patterning of trenches in 
72 Ge 28 Si layer with the overlying holes in 28 Si layer. Tak- 
ing into account that the size of all components is not 
less than 30-50 nm, one can conclude that the alignment 
could be realized by the recent progress in SEM- and 
AFM-assisted lithography. 



III. TWO-QUBIT OPERATION 

In this section, we suggest a new mechanism of entan- 
glement for distant qubits and discuss, first, the princi- 
ples of two-qubit operation. It has been shown [26] that 
two-bit gates applied to a pair of electron or nuclear spins 
are universal for the verification of all principles of quan- 
tum computation. 

Because direct overlap of wavefunctions for electrons 
localized on P donors is negligible for distant pairs, 
we propose another principle of coupling based on the 
placement of qubits at fixed positions in a quasi-one- 
dimensional Si nanowire and using the indirect interac- 
tion of 31 P nuclear spins with spins of electrons local- 
ized in the nanowire which we will call hereafter as " ID- 
electrons" . This interaction depends on the amplitude of 
the wavefunction of the " ID-electron" estimated at the 
position of the given donor nucleus ^f n (ri) and can be 
controlled by the change in the number of " ID-electrons" 
N in the wire. 

At N = 0, the interqubit coupling is totally suppressed, 
each 31 P nuclear spin interact only with its own donor 
electron. This situation is analoguous to that one sug- 
gested in the Kane proposal [5] and therefore all singlc- 
qubit operations and estimates of the decoherence time 
are valid also in our case. The difference consists in the 
method of coupling when a controlled number N of " ID- 
electrons" is injected into the nanowire. In this case, nu- 
clear spin-qubits will also interact with the spins of " ID- 
electrons" . To estimate the intensity of this interaction 
we need to calculate ^„(r i ). In the below calculation wc 
believe that the donor potential does not influence the 
distribution function \I/ n (r) of " ID-electrons" because it 
is screened by the donor electron on the relatively short 
distance of order of ae ~2.5 nm and ob is small comared 
with the wire length I. 

Let the interqubit distance be r = 200 nm, one order of 
magnitude larger than in the Kane proposal [4] . To real- 
ize the coupling between these distant qubits, we suggest 
fabricating a Si nanowire of length / = 400 nm and place 
P donors at distances r\ = (1/4)1, and — (3/4)? (see 
Fig. 3). For N = 1 and N = 3, the functions *„(r) be- 
long to the energy levels E n (n = 1, 2) because each level 
contains two electrons with opposite spin. The highest 
amplitude of ^ n (r) evaluated at the positions of the nu- 
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clear spin qubits n and r 2 is realized at N — 3 (n = 2). 
In this case, the intcrqubit coupling is maximal. 

To calculate the coefficient of the hyperfine interaction 
between nuclear and electron spins, we consider follow- 
ing [27], a system consisting of electrons confined by a 
potential V (r) and two nuclear spins. We suppose that 
the nuclear spins are located far enough from each other 
so that the direct nuclear spin interaction is negligible. 
The contact hyperfine interaction between electrons and 
nuclear spins leads to an indirect nuclear spin interac- 
tion. Let the quantum wire have finite length I in the 
x-direction with the two nuclear spins located at f\ and 
r*2 in a magnetic field H directed in the z direction. The 
Hamiltonian of the system is given by 

H = H + H x = -L (p+-A) 2 + V (r) 

8ir 

- 2.9^ B crH + —HBlnKt-iab {r-n) 

&7T 

+ -^-HBjntihvS (r - r 2 ) (1) 

where Ho is the Hamiltonian of the electron in the meso- 
scopic structure in the magnetic field, Hi = h[^ + h[ 2 ^ 
is the contact hyperfine interaction, m e is the electron 
mass, A is the magnetic-field potential, fis is the Bohr 
magneton, 7„ is the nuclear gyromagnetic ratio, 1^2 and 
a are nuclear and electron spins, and f 1;2 is radius vector 
of the nucleus. 

The effective nuclear spin interaction energy calculated 
by second-order perturbation theory is [28]: 

E= £ ^ ^ ^(l-//)+cc 

(2) 

Here, "fi and Ei are the eigenfunctions and eigenvalues 
of Ho, and fij is the electron distribution function in the 
initial and final states. We will use expression (2) to find 
the effective interaction between nuclear spins. 

We suppose that the transverse dimensions of the 
quantum wire are much smaller than its length and the 
cyclotron orbit of the electron. The confining potential 
is V(x,y,z) = V(x) - V 6(y)6(z), where V(x) = if 
^ x I, and V(x) = oo otherwise. 

The wave function should satisfy the following bound- 
ary condition: (0) = ^ n (I) = 0. The solution has the 
form 




where 5y and Sz are the transverse dimensions of the 
electron wave function. 



Let us consider the problem at T = 0. In this case the 
electron distribution function / is 1 for the filled states 
and for the empty states. Inserting this wave function 
into Eq. (2) and assuming that the Zeeman splitting 
energy in (4) is much less than the energy gap between 
levels with different n, one obtains the following expres- 
sion for the nuclear spin interaction constant A [27]: 

A / 32 MB7n fi \ 2 sin 2 [gMri/Q] sin 2 [^n(r 2 /l)] 
\ 3lSySz J 9HbH 

(5) 

It is seen, that at N = 3, A is maximal for both qubit 
positions n = (1/4)/ and r 2 = (3/4)Z. Let us estimate 
the error caused by unavoidable fluctuations in the posi- 
tions of nuclear spins in the wire. In our device, the size 
of 30 Si-spot is 30 nm. Therefore, one can expect that 
the position of NTD-introduced P donor will fluctuate 
around the central point within ±15 nm, which is about 
1/40 of the total wire length (600 nm), Ar/l < 0.025. 
In our model, the coupling is realized in the case when 
the wave function of "ID-electron" ^(r) has the maxi- 
mal value at places of the qubit location, where the space 
derivative d^f/dx is close to zero. This makes the cou- 
pling insensitive to the form of distribution function and 
to small fluctuations in the qubit positions. 

Thus, the above consideration shows that the indirect 
coupling is maximal at N = 3, while at N — 0, the 
coupling is totally suppressed. This means that in our 
model, the entanglement between two distant qubit can 
be effectively controlled by the proper variation of N. 

IV. SCALABILITY 

Scalability is the one of the most important require- 
ments of the quantum computer proposals [1]. We sug- 
gest below the schematics of the device architecture 
(Fig. 4) which satisfy the scalability requirements. It 
is worth mentioning that the above method of coupling 
opens a way to avoid the problem connected with the 
break in the one-dimensional array of qubits. This prob- 
lem is inevitable in all proposed technologies. In the 
method of coupling based on the direct overlap of electron 
wave functions [4, 6], this requirement is crucial because 
any break in the one-dimensional array of qubits stops 
the entanglement along the array and make quantum 
computation impossible. In our model, entanglement can 
exist even in the unlikely case of two or more breaks in 
the qubit array one after another, because indirect cou- 
pling can passing over wrong sites by the proper choice 
of the nanowirc length I and the number of electrons TV 
in the wire. 

Figure 4 shows the schematics of the device architec- 
ture which allows one to vary I and N. The device con- 
sists of a 28 Si nanowire with an array of 30 Si spots. Each 
spot is supplied by the overlying A-gate, the underlying 
SD-shanncl and the lateral N-gate. After NTD, P donors 
will appear in most of the spots (which transforms these 
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FIG. 4: Schematics of a 28 Si nanowire L with an array of 30 Si 
spots (qubits and non-qubits after NTD). Each spot is sup- 
plied by overlying A-gate, underlying SD-channel and lateral 
N-gate. This device architecture allows to realize an indirect 
coupling between any distant qubits (see text). 



spots into qubits) and not appear in other spots (non- 
qubits). This situation is shown schematically in Fig. 4 
where one assume that the spots 3 and 4 are non-qubits 
("0-spots") and one need to provide coupling between 
qubits 2 and 5. For this purpose, it is necessary to con- 
nect the gates N 2 , N 3 , N 4 and N 5 . The negative voltage 
applied between other N-gates and the wire contact L will 
lead to pressing-out " ID-electrons" from all correspond- 
ing areas and formation of the nanowire with I = 800 nm 
between the sites 2 and 5 only (shown in grey in Fig. 4). 
The coupling between qubits 2 and 5 will be realized via 
injection in the wire of the necessary number of electrons 
N, using the positive voltage applied to the gates N 2 -N 5 . 



In this particular example, the maximal coupling will be 
realized at N = 7, while at N = 0, the coupling will be 
totaly suppresed. 



SUMMARY 



A new method of a quantum information processing is 
suggested based on the precise placing of 31 P nuclear spin 
qubits in a quasi-one-dimcnsional 2 Si nanowire. The 
fabrication method includes the isotope engineering of Si 
and SiGe layers followed by the neutron-transmutation 
doping of the obtained structures. The interqubit cou- 
pling is based on the indirect interaction of 31 P nuclear 
spin with the spin of electron localized in the nanowire. 
The advantage of the proposed method of operation con- 
sists in its ability to control the coupling between distant 
qubits and even between qubits separated by non-qubits 
nodes in a one-dimensional array. The last feature al- 
lows one to develop the basic unit and the scalable ar- 
chitecture of a nuclear spin-based quantum computer. 
Numerical estimates show the feasibility of the proposed 
methods. 
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